An effective method of calculating the non-Markovianity J\f for single channel open 

systems 
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We propose an effective method which can simplify the optimization of the increase of the trace 
distance over all pairs of initial states in calculating the non-Markovianity TV for single channel open 
systems. For the amplitude damping channel, we can unify the results of Breuer et al. [Phys. Rev. 
Lett. 103, 210401 (2009)] in the large-detuning case and the results of Xu et al. [Phys. Rev. A 81, 
044105 (2010)] in the resonant case; furthermore, for the general off-resonant cases we can obtain 
a very tight lower bound of N. As another application of our method, we also discuss N for the 
non-Markovian depolarizing channel. 
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I. INTRODUCTION 

Inevitable interaction with external environment may 
lead to the phenomenon of decoherence for open quan- 
tum systems. In general, the assumption of Marko- 
vian approximation usually was applied to dynamical 
evolution of system. However, recently people found 
that the Markovian processes without memory and non- 
Markovian processes with memory can lead to distinctly 
different effects on decoherence and disentanglement of 
open systems. Thus the non-Markovian dynamics have 
become increasingly important and are extensive studied 
in both discrete variable [H-Q and continuous-variable 
[ilil - [l5j systems. It has been found that the non- 
Markovian effect of environment can extend significantly 
the entanglement time of the qubits [l[ and has been ex- 
perimentally observed [7j|. Because of the importance of 
the non-Markovian effect of environment, some authors 
[IB - fl9j have developed some measures to detect the non- 
Markovianity of open systems from different points of 
view. Breuer et al. [l7[ proposed a computable measure 
JV to detect the non-Markovianity of open systems. The 
idea is based on the distinguishability of quantum states 
which results from information flow between the open 
system and its environment. Rivas et al. [i~c| also pro- 
posed a measure of non-Markovianity which is based on 
the fact that negative rates are linked to whether entan- 
glement between the system and an ancilla can increase. 
Lu et al. (l9| defined a measure of non-Markovianity us- 
ing quantum fisher information flow. As we know that 
evaluation of J\f requires optimization of the total in- 
crease of the trace distance over all pairs of initial states, 
which is very difficult to accomplish. In Ref. [ijj] Breuer 
et al. considered a two- level system interacting with a 
reservoir which possesses the Lorentzian spectral prop- 
erty, and in the large detuning case they found by nu- 
merical simulation the pair of initial states pi(0) = |e) (e| 
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and P2(0) = \g)(g\ which make the optimization of the 
total increase of the trace distance (It should be noted 
that in Ref. [l7[ they use |+) and |— ) instead of |e) 
and | g) to represent the excited and ground states re- 
spectively). Later, for the same two initial states Li et 
al. [2(| obtained the analytical expression of the trace 
distance D(pi(t), P2(t)) = \h(i)\ 2 . Very recently for the 
same model in the resonant case, Xu et al. [21] found 
two different initial states which make the optimization 
of the total increase of the trace distance by using an 
analytical method and the corresponding trace distance 
D(pi(t),p 2 (t)) = \h(t)\. It should be noted that in Ref. 
[2ll ] they use b(t) instead of h(t) to represent the ampli- 
tude damping of the excited state |e) . It is intriguing that 
these two results are quite different, that is, in these two 
cases the two initial states which make the optimization 
are different, and the corresponding trace distance are 
also different. Because the optimization is very difficult 
to accomplish, until now we have not seen any reports 
about the non-Markovianity for the general off-resonant 
case. In Refs. [13, HJ the authors separately dealt with 
the optimization by different methods, one is numerical 
and another is analytical. For the same model only for 
different parameter regimes the results are quite differ- 
ent, therefore an unified understanding of these results is 
in demands. 

In this paper, we propose an effective method which 
can easily optimize the increase of the trace distance 
over all pairs of initial states in calculating the non- 
Markovianity N for single channel open systems. For the 
amplitude damping channel, we analytically re-derive the 
results of Ref. 1171 in the large- detuning case and the re- 
sults of Ref. [211 ] in the resonant case; furthermore, in 
the general off-resonant cases we can obtain a very tight 
lower bound for M. Thus an unified understanding of 
the results of Ref. [13] and Ref. [2l| is given. As an- 
other application of our method, we also discuss N for 
the non-Markovian depolarizing channel. 

The paper is organized as follows. In Sec. II, we intro- 
duce our method. The non-Markovian amplitude damp- 
ing channel and the non-Markovian depolarizing channel 
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are examined by our method in Sec. III. Finally, we give 
the conclusion of our results in Sec. IV. 



II. METHOD OF OPTIMIZING THE TOTAL 
INCREASE OF TRACE DISTANCE 

Recently, Breuer et al. [TzJ proposed a measure to de- 
tect the non-Markovian behavior of quantum processes in 
open systems based on the distinguishability of quantum 
states. The trace distance D describing the distinguisha- 
bility between the two states is defined as (22| 



D(pi,p 2 ) = 2 tr IPi - Pal 



(1) 



where \M\ = VM^M, and < D < 1. If D = 0, the 
two states are the same, and if D = 1, the two states are 
totally distinguishable. 

Considering a quantum process 3>(i), p(t) = <§>(t)p(0), 
where p(0) and p(t) denote the density operators at time 
t = and at any time t > 0, respectively, then the non- 
Markovianity M is defined as 



7V($) 



max 

Pl,2(0) 



dta(t,p 1<2 (0)), 



(2) 



cr>0 



where er(i, pi, 2 (0)) is the rate of change of the trace dis- 
tance defined as 



the trace distance for arbitrary time interval [r™ 
That is 



max [£>(pi(t, 

Pi, 2(0) 



),P2(C aX )) 

,p 2 (C in ))]. 



(5) 



Apparently it is much easier to find the two initial states 
which makes the maximization in Eq. (5) than to find 
the two initial states which makes the summation in Eq. 
(4) maximal. Then for specific non-Markovian channel 
we try to prove that the two initial states we found can 
also make the summation in Eq. (4) maximal. Generally 
speaking it is not easy to prove this. If it can not be 
proved, we still believe that the pair of initial states which 
make the increase of the trace distance in single time 
interval maximal, will also make the optimization of the 
summation in Af '. Of course this is not rigorous. If we 
are strict enough, at least in this case we can find a lower 
bound of JV($), and we argue that this lower bound is 
tight. We will show in the following that this method is 
very effective. 



III. APPLICATIONS 

Based on the above idea, we can calculate the non- 
Markovianity J\f for the non-Markovian amplitude damp- 
ing channel and the non-Markovian depolarizing channel. 



a(t, Ph2 (0)) = -D( Pl (t),p 2 (t)). 



(3) 



As we known that a(t, pi, 2 (0)) < corresponding to all 
dynamical semigroups and all time-dependent Markovian 
processes, and a process is non-Markovian if there exists 
a pair of initial states and at certain time t such that 
c(^Pi,2(0)) > 0. Physically, this means that for non- 
Markovian dynamics the distinguishability of the pair of 
states increases at certain times. 

In view of Eq. (3), the non-Markovianity M also can 
be written as the following form 

AA($) = max £ [D( Pi (t™™), p 2 {r^)) 

Pi,2(0) n (4) 

-D( Pl (r^),p 2 (r^))}, 

where r™ ax and r™ in correspond to the time points of 
the local maximum and minimum of D(pi(t), p 2 (t)), re- 
spectively. A/"($) can be calculated as follows: One first 
derives the increment of the trace distance over each time 
interval [r™ ln , r™ ax ] for any pairs of initial states, then 
sums up the total contributions of all intervals, finally 
performs the maximization for all pairs of initial states. 

Generally speaking it is very difficult to make the max- 
imization in Eq. (4). In this paper we want to find easy 
ways to maximize the increase of trace distance. Our idea 
is like this: First we find the two specific initial states 
which make the maximization of the quantity J\f n ($>) at 
each time interval. The A/" n ($) is defined as the differ- 
ence between the local maximum and local minimum of 



A. Non-Markovian amplitude damping channel 

We consider a two-level system (qubit) interacting with 
a zero temperature reservoir. The Hamiltonian of the 
total system under the rotating wave approximation is 
given by (h = 1) 



N 



fe=i 



N 

-^2(9kV-ai+g* k a + a k ) 7 (6) 

k=l 



where a + = \e)(g\ and <r_ = |g)(e|, are the Pauli raising 
and lowering operators for the two-level system, respec- 
tively, luq is the Bohr frequency of the two-level system, 
Sfc and a\ are the annihilation and creation operators for 
reservoir mode fc, is the frequency of the mode k of the 
reservoir, and g k is the coupling constant. The Hamilto- 
nian of Eq.(6) can describe various systems. For concrete 
discussion we take a two-level atom interacting with the 
reservoir formed by the quantized modes of a high-Q cav- 
ity. The dynamics of the reduced density matrix for the 
two-level atom can be written as 23] 



P 8 (t) 



p S ee^)\h(t)Y 



p s eg m(t) 

S ^il/,m2 



p S eg *(0)h*(t) 1-^(0)|^)I 



(7) 



in the basis {|e), \g}}, where the superscript S represents 
the two-level atom. Corresponding h(t) denotes the am- 
plitude of the upper level |e) of the atom initially pre- 
pared in |e) and satisfies the following integrodifferential 
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equation 



d , , . 
dt h ® 



dtif(t-ti)h(h), 



(8) 



where the kernel f(t—ti) — J duj J(u) exp[i(wo— w)(t— ti)] 
is related to the spectral density J(w) of the reservoir. 
The model describes the damping of a two-level atom in 
a cavity. In this paper we restrict ourselves to the case 
that the atom-cavity system has only one excitation, and 
suppose that J{uj) takes the Lorentzian spectral density 
23] with detuning, namely 



J(u) 



1 



7oA 2 



2tt (w - 6 - w) 2 + A 2 



(9) 



Here 6 = ujq — uj c is the detuning of the center frequency 
of the cavity w c and the Bohr frequency of the two-level 
atom wo, the parameter A defines the spectral width of 
the coupling, which is associated with the reservoir corre- 
lation time by the relation tb = A -1 and the parameter 
7o is related to the relaxation time scale tr by the rela- 
tion tr = 7 ~ 1 . Therefore the analytic expression of h(t) 
can be obtained as 

h(t) = e~ {x ~ lS)t/2 [cosh{dt/2) + (A — iS) smh{dt/2)/d] 

, (10) 

with d = y/(X - i6) 2 - 2 7o A. 

Based on the Hermiticity, and unit trace of a physical 
density matrix, any pair of initial states can be defined 
as 21 j 



pf (0) = 
P$(P) = 




(11) 



with \/3\ 2 < a(l — a), |^| 2 < p(l — p) corresponding to the 
scmipositivity of a density matrix, /3, v £ C, < a, p < 1, 
and a, p G M. Thus, the evolution of the corresponding 
density matrix can be obtained 



pf («) = 

(4(t) = 



a\h(t)\ 2 j3h(t) 

f3*h*{t) l-a\h(t)\ 2 
p\h(t)\ 2 vh{t) 

v*h*{t) l~p\h(t)\ 2 



(12) 



The combination of Eqs. (1) and (12) immediately pro- 
vides the expression of the trace distance at any time 
t > 



D(p?(t),pi(t)) = V\h(t)Ha-py + \h(tW\(3-v\i 7 

_ (13) 
which has been obtained in Ref. [21(. It is noted that 
the maximization of trace distance in the resonant case 
has been given in Ref. [21], however, their method can- 
not be extended to the general off-detuning case because 
the condition |/i(T™ m )| = can not always be guaranteed 
at each local minima. Now using Eq. (5), we can eas- 
ily achieve the maximization of Af n analytically, namely 



the optimization of the trace distance difference between 
the local maximum and local minimum in the time in- 
terval [r™ m , r™ ax ] by choosing two specific initial states. 
According to Eqs. (5) and (13), N n can be written as 



max Whir™ 

P1,2(0) U 



0M^r x )l 2 (a-M) 2 + l/3-H 5 



- IMC in )l VIM^ min )l 2 (« - m) 2 + \I3 - v\ 2 ]. 

(14) 

When t = 0, |h(Q)| = 1, and from Eq. (13) D = 
y/{a — /i) 2 + \/3 — v\ 2 < 1. Furthermore, the condition 
is equivalent to these parameterized conditions (a — p) = 
rcosfl, /3 - v = re^saxO (r < 1, Q e [0,2tt] and </> G 
[0,7r]). Substituting these parameterized conditions into 
Eq. (14) and considering that the maximization is over 
all pairs of initial states, we c an obtain that the max- 
imization condition requires a/ (a — p) 2 + jj8 — i^| 2 = 1 
corresponding to r = 1. Then the problem becomes the 
maximization of the following N°{9) 



M°(6) = \h{T™*)\^\h{T™^)\ 2 cos 2 e + sm 2 9 



\K T n m )\y \H T n in )\ 



■ + sir 



From the following equation 

dQ 



= 



(15) 



(16) 



we can obtain the extrema, which are M°i = A 2 — B 2 
when 6 = 0; Af° 2 = A 



K 3 = A 2 



l-A 2 



(B 2 -1)(A 2 +B 2 -1 

arccos 



B when 6 = tt/2 or 6 = 3?r/2; 
- B 2 



l-B 2 



(A 2 -1)(A 2 +B 2 -1) 



when 9 

A=\h(r™)\,B = \h(T, 
represented as 



-A^+2A^-A''+B^-2B^+B6 1 ,. rlnprp 
(A^-lHB^-lHAi-A^+B^-B^J' wnere 

min ^| and A > B. So Af n can be 



t{AC,A/; o 2 ,A/; 3}. 



(17) 



From numerical calculation, we find that for any A and 
B satisfying < B < A < 1, Af% 3 is always less than M^ x 
and Af° 2 , so Af n — max.{Af° 1 ,Af° 2 }. From the definitions 
of Af°x and N° 2 we can obtain that when < A + B < 1, 
K°i < when A + B = 1, Af° t = M° 2 ] when A + B > 

(i) In the resonant case, that is S = 0, it is obvi- 
ous that A + B = |/i(C ax )| + \h(T™ in )\ < 1 because 
B = |MC in )l = at r" in = 2[mr - arctan(d't/2)]/d' 
with n=l,2,3,..., and d' = \/\\ 2 — 27oA|. So in any time 
interval [r™", t™% N n - nua^,^} = N° 2 = 
A- B = |/i(C ax )|. For Af° 2 , 9 = tt/2 or 3tt/2, and the 
two initial states correspond to a = p and |/3 — v\ = 1. 
Because in this case Af n reaches its maximum for any n, 
the same two initial states, that is a — p and \0 — v\ = 1, 
are also the two initial states which make the summation 
in Eq. (4) maximal. It is worth noting that these con- 
ditions a = fi, \j3 — v\ = 1 together with \f3\ 2 < a(l — a) 
and \v\ 
a = p = 1/2 
tained in Ref. 



< p(l — p) are equivalent to the conditions 
= \v\ = 1/2 and \fi - v\ = 1 ob- 
which can be explained as follows. 
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Our conditions can be changed into a — fx, \/3 — v\ = 1, 
(a-l/2) 2 + |/3| 2 < (1/2) 2 and {jjl- 1/2) 2 + \v\ 2 < (1/2) 2 . 
From the geometric point of view, the new conditions in- 
dicate that the two points (a, \/3\) and (/x, \u\) is in (or on 
the circumference of ) the same circle 

0r-l/2) 2 + M 2 = (l/2) 2 (18) 

with x G M and y G C. It is easy to check that the 
conditions a = [i, and \/3 — v\ = 1 are just a = /i = 1/2, 
|/3| = |^| = 1/2, |/3 — ^| = 1. In summary, in the resonant 
case our results are consistent with the results of Ref. 
21]. 

(ii) In the off- resonant case, that is 5 ^ 0, A + B = 
|ft(r™)| + |MT™ in )| may be less than 1, equal to 1, or 
more than 1 depending on the values of 70, A and 5. We 
have proved that depending on the value of A + B there 
are only two pairs of initial states which make the maxi- 
mization of Af n for each n. Next we give our effective and 
practical method to calculate Af for any fixed parameters 
70, A and 8. From Eqs. (2) and (3) we can use the two 
pairs of initial states to obtain Af\ and A2 respectively, 

= L>o dt °^ 6 = °)> ^ = L 2 >o dt(J ^ 9 = */ 2 OT 
3ir/2). Correspondingly, the expressions of a\(t,8 = 0) 
and <T2(t,6 = ir/2 or 3ir/2) are given by 

<Jx(t) = e~ At {/i[cosh(crf) — cos(6t)]+^sinh(at) — £ sin(fet)}, 

(19) 

e~ 2 £ {^[cosh(ai) — cos(bt)} + ^sinh(ai) — £ sin(Zrf)} 

&2{t) = , 

2 -y/ 77 sinh(ai) + xsin(6i) + Kcosh(ai) + <;cos(bt) 

(20) 

where a and b denote the real part and imaginary part of 
d respectively; \x = ( Aa 2 — A6 2 — X5 2 — A 3 — 2abS) , v — 

1 ^{a 3 +ab 2 -\ 2 a+a5 2 + 2b5\), £ = 1 ^ I (b 3 +ba 2 + \ 2 b- 

bS 2 + 2aSX), n = 2^p( 2aA - 2M )> X = 5Rp(2&A + 2aS )> 
K = 2^r(A 2 + (5 2 + a 2 + 6 2 ), ? = ^( a 2 + b 2 - X 2 -S 2 ); \d\ 
denotes the absolute value of d. It is worth noting that 
Eq. (19) has been obtained in Ref. [20|. We can not 
prove but we believe that one of the two pairs of initial 
states we found can also make the optimization in the 
summation of Eq. (4), thus Af — max{A/i , A2 } ■ If we 
are strict enough, at least it is a very tight lower bound 
(TLD) for the genuine Af, 

A/"tld = max{A/"i , A/2}. (21) 

We plot A/"tld, Ai, and A/2 as functions of detuning 
S for A = O.I70 in Fig. 1, and A/"tld, A/i, and A/2 as 
functions of A for S = O.I70 in Fig. 2. From Fig. 1 we 
can see that there exists a critical point for 6 at which the 
pair of initial states change from 8 = ir/2 or 8 = 3ir/2 
to 8 = 0. More specifically when 5 < 5 C , A/"tld = A/2 
corresponding to the two initial states 8 — it/ 2 or 8 = 
37r/2; when 5 = 8 C , Atld = A/i = A2 corresponding to 
the two initial states = 0, and 6> = n/2 or 37r/2; when 
5 > 5 e , A/"tld = A/i corresponding to the two initial 




FIG. 1: (Color online) Atld, A/i and A/2 as a function of (5, 
A = 0.l7o. 




FIG. 2: (Color online) A/tld, A/i and A/2 as a function of A, 
S = O.I70. 



states 8 — 0. Similarly, it can be seen from Fig. 2 that 
there also exists a critical point for A at which the pair of 
initial states change from 8 = to 8 = tt/2 or 8 = 3ir/2. 

It is worth noting that in the large detuning case by 
numerical simulation Breuer et al. [171 ] have performed 
the optimization of the total increase of the trace dis- 
tance. In the large detuning case from Eq. (10) we know 
A + B = |/i(-C ax )| + IMC in )l > 1 for any n = 1,2,3,.... 
Therefore, we can obtain M n = max{A/"° 1 , A/"° 2 } = Af°± = 
A 2 -B 2 = |/i(C ax )| 2 - IMC in )| 2 for any n and the two 
initial states satisfy a — fi = 1 and |/J — v\ — corre- 
sponding to 8 = 0. Because the two initial states (8 = 0) 
make the optimization of Af n for any n, they also make 
the optimization of the summation in Eq. (4) . Similar to 
the resonant case, the optimization conditions including 
a-fi = 1, \p-v\ = 0, \Ci\ 2 < a(l-a), and \v\ 2 < ti(l-fi) 
can be simplified into a = 1, /J, = and /3 = v = 0. 
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Evidently, the two initial states are pi(0) = \e)(e\ and 
P2(0) = which has been obtained in Ref. 17]. It is 

noted that the condition is obtained by numerical simula- 
tion in [17], but here we obtain the same condition using 
an analytic method. Then we can also obtain the trace 
distance D = |fe(t)| 2 for these two initial states which is 
given by Ref. 20]. In summary, in the large-detuning 
case our results are consistent with the results reported 
inRefs. SUS- 



B. Non-Markovian depolarizing channel 

As another application of our method, we consider the 
non-Markovianity Af for a non-Markovian depolarizing 
channel. The dynamical property of this system and in 
particular the conditions of complete positivity of the 
map corresponding to a master equation have been stud- 
ied in detail by Daffer et al. in Ref. [24j. For this model 
the time-dependent Hamiltonian, that corresponds to a 
two-level system subjected to random telegraphic noise, 
is 



H(t) = hJ2Ti(t)ai, 



(22) 



where I\(i) = aiUi(t) are independent random variables, 
and <Ji are the usual Pauli operators. rii(t) has a Pois- 
son distribution with a mean equal to i/2ri, while ai 
is an independent coin-flip random variable taking the 
values ±a;. For the time-dependent Hamiltonian of Eq. 
(22), the corresponding equation of motion for the den- 
sity operator is govern by the von Neumann equation 
p = —i/h[H(t),p] = — iY^k^k{t)[(Jk, p] which has the fol- 
lowing formal solution 



p(t) = p(0) - i 



JO u 



(s)[(x k ,p(s)}ds. (23) 



Substituting the formal solution Eq. (23) into the von 
Neumann equation and performing a stochastic average, 
one can obtain the following memory kernel master equa- 
tion 



m = 



y^ e -{t-t')/r k 



al[a k ,[a k ,p(t')]}dt', (24) 



where the correlation function of the random telegraph 
signal (r j (t)r fc (t')) = ale-^-t'^^Sjk contributes to the 
memory kernel. It has been pointed out [HJ that the sys- 
tem density operator with an exponential memory kernel 
obeys a homogeneous Volterra equation after averaging 
over the reservoir variables, and also proved that when 
two of at are zero, only one direction having the noise, 
the map &(p) can be written as Kraus operator form [25| . 
namely 



o'{t) = <S> t {p) = Y J A\pA k 



(25) 



For simplicity, in this paper we only consider the case the 
z direction with noise, x and y directions without noise. 
Therefore, A 1 = 0, A 2 = 0, A 3 = y/[l - A(v)]/2a 3 , and 
At = + A{v)]/2I, where A(v) = exp(-i/)[cos(/xz/) + 
&m(pv)/p\ with p — y/ (4<2t) 2 — 1, and v = t/2r is a di- 
mensionless time, a is the coupling strength of the system 
with the external environment while r determines which 
frequencies the system prefers most. For convenience we 
let A = 1/r, then A{v) can be rewritten as 



A(t) = 

' e- At / 2 [cosh(f ) + j sinh(f )] (16a 2 < A 2 ) 



e~ At / 2 [l + f] 
t e- A */ 2 [cos(f ) + f sin(f )] 



(16a 2 = A 2 ) 
(16a 2 > A 2 ), 



(26) 



where e — y/\16a 2 ~ A 2 |. 

Using the same two initial states of Eq. (11), from 
Eq. (25) we can obtain the evolutions of the two density 
matrixes, respectively 



pf(t) 

pi(t) 



a /3A(t) \ 
P*A*(t) I- a) 
p vA(t) \ 



v*A*(t) 1 - p J • 
Therefore the trace distance can be obtained 



D(p?(t),p§(t)) = y/(a-tf + \A(t)\*\l3-v\*. 
From Eqs. (5) and (28), Af n can be expressed as 



Af n = max [y/ (a - p) 2 

Pl,2(0) 



\A(Tr x )\ 2 w-v\- 



^{a-pf + \A{r^W- V \ 



(27) 



(28) 



(29) 



As in the subsection A, after parameterizing (a — p) and 
|/3 — u\, the problem becomes the maximization of the 
following N° {6) 



N°{0) = Jcos 2 e+ |A(r^ ax )| 2 sin 2 6» 



cos 2 6»+ \A(t£ 



(30) 



Then we can obtain the extrema of N°: Af° ± = when 
9 = 0; N° 2 = A - B when 9 = vr/2 or 3tt/2, where 
A = \A(t™*)1 B = |A(t™)| and A > B. Obviously, 
N n = msx.{M° 1 ,M° 2 } = N° 2 = A - B. Because the 
pair of initial states corresponding to 8 = tt/2 or 37r/2 
makes the increase of trace distance Af° maximal for any 
ri, it also makes the summation in Eq. (4) maximal. 
The condition 9 = tt/2 or 37r/2 means a — p = and 
|/3 — v\ = 1, and similar to the discussion in subsection A 
the condition together with |/3| 2 < a(l — a) and |^| 2 < 
p(l — p) is also equivalent to a = p = 1/2, |/3| = \u\ — 
1/2, |/3 — ^| = 1. The trace distance after choosing the 
two initial states can be given as 



fe=i 



D = |A(i)|. 



(31) 
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and £,p(R, t) are positive, monotonically decreasing func- 
tions under the condition of positivity of the dynamical 
maps [13] ■ Thus the quantum processes corresponding to 
Eqs.(4) and (10) of Ref. [26[ are Markovian. However, 
the model we used can also be described by the mem- 
ory kernel master equation Eq. (24), but we clearly see 
that the function A(t) which plays a central role in our 
model, isn't a monotonous function with time but is a 
damped oscillating function in some parameter regimes, 
which can be seen from Eq. (26). And in this case the 
quantum process we consider is non-Markovian. 



FIG. 3: The non-Markovianity Af as a function of i 



IV. CONCLUSIONS 



It is noted that there are some similarity between the 
non-Markovian depolarizing channel and the amplitude 
damping channel. In the non-Markovian depolarizing 
channel the maximal trace distance is a function of the 
decoherence factor |A(i)|, while in the amplitude damp- 
ing channel it is a function of the amplitude damping 
factor |/i(t)|. It is very interesting that the trace distance 
of the non-Markovian depolarizing channel for the spe- 
cific pair of initial states which make the optimaization, 
is very similar to that of the amplitude damping channel 
in the resonant case, that is for the former the trace dis- 
tance for the specific pair of initial states is equal to the 
decoherence factor |A(t)|, while for the latter it is equal 
to the amplitude damping factor furthermore, in 

both cases the two initial states which makes the opti- 
mization are the same. 

We can obtain the rate of the change of the trace dis- 
tance using Eq. (3) 



<K*,Pl,2(0)) 



|A(t)| ■ 



(32) 



Having these preparations, we can easily calculate the 
non-Markovianity Af from Eq. (2). We plot Af as a func- 
tion of a in Fig. 3. From Fig. 3 we can see that there is a 
threshold a c = A/4. When a < a c , Af = 0, which means 
that the process is Markovian. While when a > a c , Af 
increases with a, which means that the process is non- 
Markovian. This can be easily understood: because a 
represents the coupling strength of the system and the 
reservoir, it is obvious that the non-Markovianity Af will 
become more and more strong with the increasing of a 
in the non-Markovian re gim e. It is worth noting that 
recently Mazzola et al. [26| elucidated that the mem- 
ory kernel master equation does not ensure the presence 
of non-Markovian behavior in the evolution of dynam- 
ics. For example, in Ref. [26[ they have shown that the 
non-Markovian behavior does not appear in the mem- 
ory kernel master equations Eqs.(4) and (10) of the same 
reference. £m(-R, i) an d £p(R,i) are two functions cor- 
responding to two models Eq. (4) and Eq. (10) of Ref. 
26] respectively, which play the central role in the dy- 
namics of the system. It has been proved that £m (R, t) 



As we know that the definition of the non- 
Markovianity Af needs an optimization over all pairs of 
initial states, and generally it is very hard to do this. In 
this paper we proposed a method which can simplify this 
optimization. The main idea is like this: First we find 
the pair of initial states which make the maximization 
of the difference between the local maximum and local 
minimum of the trace distance in arbitrary nth time in- 
terval, then we try to prove that this pair of initial states 
can also make the optimization of the summation over all 
pairs of initial states in calculating the non-Markovianity 
Af. Using this method we have analytically obtained the 
pair of initial states which make the optimization in Eq. 
(4) and the corresponding trace distance for the ampli- 
tude damping channel in both the resonant case and the 
large-detuning case, and unified the results of Breuer et 
al. [TtJ and Xu et al [2l|; and we also have analytically 
obtained the pair of initial states which make the opti- 
mization in Eq. (4) and the corresponding trace distance 
for the non-Markovian depolarizing channel. Generally 
speaking it can not always be proved that the pair of 
initial states which make the maximization of the nth 
difference between the local maximum and local mini- 
mum of the trace distance, also make the optimization 
of the summation in Af. But here we argue that the pair 
of initial states which make the increase of the trace dis- 
tance maximal in one time interval, will also makes the 
optimization of the summation in Af. We can not prove 
this in this paper, and this needs further investigations. 
If we are strict enough, at least we can obtain a very 
tight lower bound for Af. For example, for the amplitude 
damping channel, we have found only two pairs of initial 
states depending on the system parameters which make 
the increase of the trace distance maximal in arbitrary 
single time interval. Then it is easy to calculate Af for 
these two pairs of initial states, and the larger one is the 
tight lower bound of Af. In this paper we mainly focus on 
the single channel case and the generalization to the gen- 
eral multi-channel case may be a more challenging task. 
In one word we have simplified the problem of finding a 
pair of initial states which optimizing the summation of 
the increase of trace distance in many time interval into 
a problem of optimizing the increase of trace distance in 



7 



just one single time interval. 



and 11075013). 
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